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ABSTRACT. By modifying Beukers' proof of Apery's theorem that ^(3) 
^ is irrational, we derive criteria for irrationality of Euler' s constant, y . For 

^ n > , we define a double integral and a positive integer 5„ , and prove 

CN that with = LCM ( 1 , . . . ,n) the following are equivalent. 

1. The fractional part of log5„ is given by {\ogSJ = /„ for some n . 

O 2. The formula holds for all sufficiently large n . 

3. Euler's constant is a rational number. 

A corollary is that if {log5„} > 2~" infinitely often, then y is irrational. 

^ Indeed, if the inequality holds for a given n (we present numerical 

^ evidence for \ <n< 2500) and 7 is rational, then its denominator does 

^ not divide ii2n(^")- We prove a new combinatorial identity in order to 

^ show that a certain linear form in logarithms is in fact logS,,. A by- 

CN product is a rapidly converging asymptotic formula for y , used by 

^ P. Sebah to compute 7 correct to 18063 decimals. 

• 
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1. INTRODUCTION 



Since R. Apery's startling 1979 proof [1], [7] that ^(3) is irrational, there have 
been attempts to extend it "up" to ^(n) = Lr=i^" n>3, with recent progress by 

K. Ball and T. Rivoal [4] and W. Zudilin [13]. In the present paper, we go "down" to 
Euler's constant. 



where Hf^ = YJl=yk ' is the A'^* harmonic number, and find criteria for it to be 
(ir)rational. Note that y can be thought of as "^(1)," and as such sits on the fence 
between rational values of ^(n) for n < 1, and irrational values (presumably) for n > 1. 
Defining the double integral 



which is inspired by F. Beukers' celebrated integrals for ^(2) and ^(3) in his elegant 
1979 proof [2] of Apery's theorem, and denoting by 5„ the integer product 



where = LCM ( 1 , . . , we prove the following necessary and sufficient conditions 
for rationality of y . (Of course, their negations are then criteria for irrationality of y .) 

Rationality Criteria for y . The following are equivalent: 

(a) The fractional part of logS^ is given by {logS^} = I„ for some n . 

(b) The formula holds for all sufficiently large n. 

(c) Euler's constant is a rational number. 



(1) 



7 := \im(Hj^-\ogN), 



(2) 




(3) 




k=l i=0 j=i+l 
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We first discovered the Rationality Criteria by modifying the rational function 
and infinite series in Y. Nesterenko's 1996 proof [6] of Apery's theorem. In this 
approach, we defined to be the sum of the series 



V=n+l 



y \x(x + l)---(x + n) 



dx 



rather than the double integral (2). Remarkably, the two definitions agree. Although we 
do not know a change of variables transforming one into the other, it turns out that 
evaluating the integral yields the same result as summing the series, namely, 

(5) i„ = (";:)r+L„-A„ 

where 

(6) L„:=d,-^\ogS„ and A„:=£(j'/f, 



1=0 



n+i 



(Note that does not actually involve d2n, since it cancels out from (3).) We prove (5) 
for (2) in §2, but defer the proof for (4) to [1 1], since this series representation of will 
not be used here. (It might, however, be useful in deciding the arithmetic nature of y from 
the Rationality Criteria.) 

After reading D. Huylebrouck' s 2001 survey [5] of multiple integrals in 
irrationality proofs, we found (2) and rederived the Rationality Criteria from it. It turns 
out that not only is the integral method more elegant, but also it yields better sufficient 
conditions for irrationality of gamma - see (ii), (Hi) in Corollary 7 - than the series 
method does, since the integral is easier to estimate. On the other hand, to prove (5) using 
the integral definition of /„ requires a combinatorial identity that is not needed if we use 

the series definition. 

An expression like that in (2) gives a new double integral for Euler's constant 



' JoJo n - 



(1 — x) dxdy 



(1 - xy) logxy 



The proof (see [10]) is similar to that of (5) and is omitted since we will not use the 
formula. 
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We prove the Rationality Criteria in §3, and in §4 we give as corollaries several 
sufficient conditions for irrationality of j ; they involve 5„ , but not . Here is the most 

stringent one. 

If {log5„} > 2~" infinitely often, then y is irrational. 
Indeed, the following stronger, finite version holds: 

For fixed n, if {log5„} > 2~" and y is <^ rational number with denominator q, then q is 
not a divisor of d2n j ; particular, \q\>2n. 



For example, since {logSjoooo} - 0.73... (according to P. Sebah [9]) and 0.73 > 2 
if y is rational, its denominator does not divide ^20000 (10000) • 



-10000 



Table 1 exhibits S„ in the form 5„ = s^''" , where ■.= djd2„ and s^^ := S^"'^''"^ 



= exp (- L,^ ) are integers ; like powers of n + ^ and 2n - /c + 1 are combined. 
TABLE 1. 



n 


Sn 


{log5„} 


1 




0.7725... 


2 


(i2^y^ 


0.4566... 


3 


(2411.538)20 


0.3446... 


4 


(4025.42185^40 


0.7212... 


5 




0.9645... 


6 




0.5546... 



Computations [9] show that {log5„}>2~" for l<n<2500. Also, they suggest that 
{log^,,} is dense in the interval (0,1) and that the cumulative average <{log5(n)}> 
= n~'£^'^j{log5^} tends to 1/2 (see Figure 1, courtesy of P. Sebah). To prove y 
irrational, though, it would suffice just to show that {log^,,} does not tend to zero. 
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HGURE 1. 

□ 75- 
0.7- 
0.65- 

<{logS(n)}> 




^■^5 500 !□□□ 15^^ 2^^^ 25^^ 

n 

In §5, we prove the asymptotic formula 

Table 2 shows the first five approximations it provides. 
TABLE 2. 



r 


= 0.577216664901.. 




n 






1 


0.5(56852819440.. 


.) 


2 


0.57(6991219550.. 


.) 


3 


0.57721(2786561.. 


.) 


4 


0.5772166(25800.. 


.) 


5 


0.577216664(353.. 


.) 



Using the formula with « = 10000, P. Sebah [9] has computed y correct to 18063 
decimals, which is precisely the accuracy 2~^''°°*'/l00 predicted. 
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Finally, in the AppendiK we establish the combinatorial identity 



k-l n 



i=0 j=k 




where Hq=0 . Then we use it to prove Lemma 2, which implies that times a certain 

linear form in logarithms (viz. (11)), a product which is a priori the log of a rational 
number, is in fact the log of an integer, namely, logS^. (A reader who wishes to skip the 
Appendix may take (11) as the definition of L„, and substitute d2„L^ for logS^ in the 

Criteria and the corollaries, although this makes their statements less elegant.) 

In a paper in preparation [12], we show that conditions stronger than (/) and (//) in 
Corollary 7 imply upper bounds on irrationality measures for y . 

I thank the referee for suggestions on exposition, and Patrick Gallagher for 
lectures on [6] and discussions of Laplace's method. 



2. EVALUATION OF THE DOUBLE INTEGRAL 



Theorem 1. The following relation is valid for all n>0: 



(7) 




where L„ is the linear form in logarithms 



(8) 




and d2„A^ e Z. 
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Proof. Fix n > . In the integrand of , we expand 1/(1 - xy) in a geometric series with 
remainder 



ifi {x{\-x)y{\-y)r 



logxy 



\k=\ ^ 



{xy) 



■xy 



dxdy 



and set 



(9) 



R 



N 



JOJ( 



in {x{\-x)y{\-y)nxy)'' 



(l-xy) logxy 



dxdy . 



In I^-Rf^, we substitute 

1 poo 

(10) =\{xyydt, 

use the Binomial Theorem and reverse the order of integration (which can be justified by 

rl rl-£ roo rl/e 

replacing by , and by , and letting e tend to 0), obtaining 

Jo Je Jo Je 



poo " ..(n\(n\^ 1 

Jo j:'o I * AiJrt {n + i + k + t)(n + j + k + t) 



0<i<j<n i=j=Oj 



^ (n + / + A: + t)in + j + k + t) 



dt. 
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In this last equality, we have used the symmetry of the expression in / and j . It is easily 
seen that for < / < j <n and N >0, after expansion in partial fractions, the inner sum 
on k telescopes (cancels) to give 



1 ^ 



1 



1 



-^'ii^ ' 

j -if^^\n + i + k + t N + n + i + k + i 



n + i + k + t n + j + k + t ^ 
Integrating these terms, as well as the squared ones for / = j, yields 



V 0<i<j<n J ^ 



(-ly+YnYnAf^' N + n + i + k 

> log 

yj jT^, n + i + k 



i^oy y o<i<j<n J yyyjyk^i 

Since n + i< In, n + i + k<n + j <2n and n is fixed, we have 

Hff+„+i =Hj^ + o(l), log(N + n + i + k) = logN + o(l) 
as N ^ oo. Therefore 



where 



(11) 



(=0 



In-RN = U']H^-2 (-l)'"ni" logA^ + L„-A„+o(l), 

• " 0<i<7<n 



0<i<i<nk=\ J * 



\og{n + i + k) 



and A„ is given in (6). Hence ^ ^ • 
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Using the identities 

2 1 (-i)"'"1:'){:)=i(:)"=f:) 

0<i<j<n i=0 

(to prove the first, put x = -l in ((1 + x) ") ^ ; for the second, which is a special case of the 
Chu-Vandermonde identity, compute the coefficient of x" in two ways), we obtain that 

as A/^^ oo 

I„-Rj,= (';) (H^-logN) + L„-A„+o(l). 

In view of (1), it only remains to show that Rj^ = <?(!), and that expressions (8) and (11) 
for L„ agree. 

For any n > 0, the quantity ^""^^"i^^^y"^^^" in (9) is bounded by 1 on [0,1]^ Using 
(10), it follows that 

SO that 0<R}^ <(N + 1)"^ — > . The following lemma completes the proof of Theorem 1 . • 
Lemma 2. The formulas for in (6), (8) and (11) all agree. 



The proof is given in the Appendix. 
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3. RATIONALITY CRITERIA FOR 7 

Lemma 3. The inequalities 1 < < 8" and < I^< 16~" hold for all n>0. 

Proof. It is known (see [8], Theorem 13) that log rf„ < 1.03883 n for all n. As 

2xl.03883<31og2, wehave 1< ^2,, <8". 

From (2), observe that /„ >0, and that t(l-t)<l/4 for implies that 

/„/16. Since /j =27 + 2 log2- 5/2 < 1/16, the lemma is proved. • 

Theorem 4. For n fixed, {\ogS^} = (ijn ift^^d only if y = p/q for some p,q sZ, 



where q 




Proof. Multiply (7) by d2„ and use (6) to write the result as 

\0gS„ - d2„ I„ = d^n K - ^2n(^„j 7 , 

where t/jw ^ ^ • Note that \ogS^ > , and that d2„ /„ e (0,1) , by Lemma 3. Since the 
fractional part of x > is the point {x} s [0,1) such that x - {x} € Z, the theorem follows. • 

Theorem 5 (Rationality Criteria for y ). The following are equivalent: 

(a) The fractional part of log5„ is given by {logS^} = Kfaf some n. 

(b) The formula holds for all sufficiently large n. 

(c) Euler's constant is a rational number. 

Proof. Since q I d2n for all n>qll. Theorem 4 implies that (c) {b) (a) (c). • 
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4. SUFFICIENT CONDITIONS FOR IRRATIONALITY OF 7 

Corollary 6. If {logS'^} > 2~" infinitely often, then y is irrational. In fact, for a given n 
the inequality implies that if p,qGX and q I tijn (^") > Y ^ V 1^- 

Proof. Lemma 3 and Theorem 4 imply the second assertion, which implies the first. • 

Corollary 7. Euler's constant is irrational if at least one of the following is true: 
(0 lim sup t" { log 5'„ } > for some positive number t <(4/e)^, 

(«) {logSJ^(e/4f"^'^"^'^\ 

lim {log5„} ^ . 

(ill) n^o. 61og2 

Proof Using = (by the Prime Number Theorem), we see that condition (/) 

(n) =^ (Hi). To show that (Hi) y ^Q, we make the change of variables u = 2x-l, 
v = 2y-l in (2), and find that 

/„=4-f c^^^^^n^^^^dudv, 

J-iJ-i(4-/(„,v))log(4/(M,v)-i) 

where f(u,v) :- (1 + u)(l + v). Since = exp(-w^ + O (w"^)) as w — > 0, it follows by 

Laplace's method (see [3], p. 322) that as n ^ 0° 



(12) i^^4-^"l±^, 

31og4 

which contradicts (Hi) if 7 e Q, since then {\ogSJ = d2„ I„ for n large. This completes 
the proof. • 

Corollary S. If {logS'^j < 16{logS„+i} infinitely often when 6?2n ~ <^2n+2' 7 
irrational. Indeed, if the inequality holds for a given n and rijn ~ ^in+i' 7 
rational number with denominator <2n + 2. 

Proof. The first assertion follows from the second, which in turn follows from 
Theorem 4, since > 16/„+i. • 
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5. AN ASYMPTOTIC FORMULA FOR 7 
Theorem 9. The following asymptotic formula holds: 



' 

Proof. This follows from Theorem 1 together with formula (12) and the estimate 
(^J-J = 0(4~"n^'^) from Stirling's formula. • 

APPENDIX. A COMBINATORIAL IDENTITY AND PROOF OF LEMMA 2 
Proposition 10. For \ <k<n, we have 



k-\ n 



(13) 



i=Q j=k 



1 



k-\ 



We require a lemma. 
Lemma 11. We have 



(14)„ 



and 



(15)„ 



k-\ 



k+i-\ 



! = 



13 



Proof. We prove (14)„ by induction on n. Let L(n,k) denote the sum on the left-hand 
side of (14)„. For n = l, we have ^ = and L(l,0) = l = H^, as required. Now take 
n > 1 . For A: = , we use the recursion 



(16) 



to break L(n,0) into two sums, in the second of which we substitute 





(n-V\ 




fn 


-n 






+ 






J. 


V j y 




J 



(17) 
obtaining 



li^ 



L(n,0) = Lin- 1,0) + -Y,(-iy~' 



7=1 



= L(n -1,0) + 



1 



Hence by induction L(n,0) = , which proves (14)„ in the case k = . 

For A:>0, we again use (16) to break L{n,k) into two sums, but now in the 
second one we set / = 7-1 and write 



L{n,k) = L{n-\,k)+ X^C-l/^'f " ./ 

= L{n-\,k) + L{n-\,k-\). 
Assuming inductively that (14)„.i holds, we have 
L{n,k) 



1 



r-{k-\) 



(n-\\ , (n-\\ , 



Kkj 



yk-\j 



{H„.,-H,) + 

I 

and an application of (17) with j replaced by k yields (14)„. 
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To derive (15)„, we start with (14)„ , set j = n-i, and replace k by n-k, and 

(:'J by (:) 

, for m = / and m = k. The result is (15)„ and the lemma is proved. • 

Proof of Proposition 10. We fix n and induct on k . Let L(k) and R(k) denote the left- 
and right-hand sides of (13). According to (14)„ with k = 0,'we have ^(1) = H^ = R(l) . 
Now observe that 



L(k + l)-L(k)= (-1)'^'"' 

j=k+i 



ykj 



1 



ykj 



1 



k — i 



for 1 < ^ < n - 1 . Applying (14)„ and (15)„, we obtain 



L{k + \)-L{k) 
which by induction proves (13). 



ykj 



{H„_,-H,) = R{k + \)-R{k). 



We recall Lemma 2 and give the proof. 
Lemma 2. The formulas for in (6), (8) and (11) all agree. 



Proof. We may write the right side of (1 1) as 



(18) 

and that of (8) as 
(19) 



n k—\ n 

2EIE 

k=\i=Q j=k 



(-1) 



i+j-l 



-\og{n + k) 



n k-lfyj^'^ 



k=\i=Q 



2^1 • (Hn-rtl,)\og{n + k). 

1 1 .• n\ I I 



By Proposition 11, the coefficient of log(n + ^) in (18) is equal to that in (19), for 
\<k<n. Hence, the four expressions for L„ in (8), (11), (18) and (19) are all equal. 
Since those in (6) and (8) evidently agree, we are done. • 
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